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Many-body-QED perturbation theory: 
Connection to the Bethe-Salpeter equation 
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Abstract: The connection between many-body theory (MBPT) — in perturbative and non- 
perturbative form — and quantum-electrodynamics (QED) is reviewed for systems of two 
fermions in an external field. The treatment is mainly based upon the recently developed 
covariant-evolution-operator method for QED calculations [Lindgren et al, Phys. Rep. 389, 
161 (2004)], which has a structure quite akin to that of many-body perturbation theory. 
At the same time this procedure is closely connected to the S'-matrix and the Green's- 
function formalisms and can therefore serve as a bridge between various approaches. 
It is demonstrated that the MBPT-QED scheme, when carried to all orders, leads to a 
Schrodinger-like equation, equivalent to the Bethe-Salpeter (BS) equation. A Bloch equation 
in commutator form that can be used for an "extended" or quasi-degenerate model space 
is derived. It has the same relation to the BS equation as has the standard Bloch equation 
to the ordinary Schrodinger equation and can be used to generate a perturbation expansion 
compatible with the BS equation also for a quasi-degenerate model space. 
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1. Introduction 
1.1. General 

What is known as the Bethe-Salpeter (BS) equation represents the complete solution of the relativis- 
tic two-body problem with important applications in various branches of physics. The equation was 
first derived by Bethe and Salpeter in 1951 [1], using the relativistic S'-matrix formalism and the anal- 
ogy with Feynman graphs, and at about the same time by Gell-Mann and Low [2], using a rigorous 
field-theoretical approach based on Green's functions. A closely related equation was discussed by 
Schwinger in his Harvard lectures already in the late 1940's [3, 4, 5, 6]. 

In interpreting the solutions of the BS equation, several serious problems were encountered, as 
discussed early by Dyson [7], Wick [8] and Goldstein [9]. Dyson was particularly concerned about 
the meaning of the wave function in relativistic quantum mechanics, a subject "full of obscurities and 
unsolved problems ". Solving the BS equation leads to a 4-dimensional wave function — with individual 
times for the two particles. This function is manifestly relativistically covariant but not in accordance 
with the standard quantum-mechanical picture. That leads to "spurious" or "abnormal" solutions with- 
out physical significance and with no nonrelativistic counterpart [10]. Another fundamental problem 
is that the BS equation does not reduce to the correct "one-body limit", when one of the particles be- 
comes infinitely heavy, as discussed by Gross and others [11, 12]. Problems of these kinds are most 
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pronounced in the scattering of strongly interacting particles but less so for bound-state systems in 
weak-coupling [13, 14, 15, 16, 17] (see ref. [6] for a review). 

The earliest applications of the BS equation appeared in atomic physics and concerned the proton 
recoil contribution to the hydrogen fine structure by Salpeter [18] and the positronium energy level 
structure by Karplus and Klein [4] . 

An important goal for the equation has been the study of strongly interacting particles, which is a 
fundamental problem in elementary-particle physics. In recent years there have been numerous appli- 
cations in QCD, dealing mainly with the quark-quark, quark-antiquark interactions, quark confinement 
and related problems [14, 19, 20, 21]. Here, the problems mentioned above are more serious, as recently 
summarized by Namyslowski [6]. 

There have also been many applications in surface and solid-state physics, ranging from electron- 
hole interactions in ion crystals [22] and studies of the two-dimensional Hubbard model [23] and 
Cooper pairs [24] to quantum dots [25]. 

The BS equation has also been applied to three or more particles [26, 27, 28], although serious 
problems have been encountered for more than three particles [29]. 

Various approximation schemes for treating the BS equation have been developed over the time. 
The simplest approximation is the "ladder approximation ", where all intermediate states evolve only in 
the forward (positive) time direction. This is a useful starting point in the strong-coupling case, where 
the standard perturbative or self-consistent approach may not converge, and this approximation is, for 
instance, the basis for the Brueckner theory of nuclear matter [30, 31, Sect. 41]. Another approach 
is the "quasi-potential approximation ", which implies that the equation is reduced to an equivalent 
3-dimensional Schrodinger equation, which can be done without loosing any rigor [14, 32]. Early 
numerical calculations in the this regime were done particularly by Schwartz and Zemach [33] and 
Kaufmann [34]. 

In atomic physics the BS equation has been applied mainly in treating positronium [35, 36] and to 
heliumlike ions, and we shall be particular concerned with the latter here. This is strictly speaking a 
three-body problem but can to a good approximation be treated — with the first Born approximation — 
as a two-body problem with an external potential. The application to heliumlike systems was pioneered 
by Sucher [37, 38] and Araki [39] in the late 1950's for deriving the leading relativistic and QED energy 
corrections beyond the Breit interaction. Later these works have been extended — largely along the lines 
of Sucher — by Douglas and Kroll in the 1970's [40] and more recently by Zhang and Drake [41, 42, 
43,44]. 

The technique developed by Drake and coworkers is presently the most accurate available in dealing 
with heliumlike systems. The wave functions used are very accurate functions of Hylleraas type, and 
the QED corrections are evaluated by means of analytical expressions up to order a 5 Ry (atomic units, 
or ma 7 in relativistic units), derived from the BS equation. The wave functions used by Drake et al. are 
nonrelativistic but certain relativistic effects are treated to all order in the "unified model" [45, 46]. The 
analysis of the BS equation are in these works based upon the Brillouin-Wigner perturbation theory 
(BWPT). 

A different and in some aspects more versatile approach to the many-body problem is the proce- 
dure known as the many-body perturbation theory (MBPT). This is based upon Rayleigh-Schrodinger 
perturbation theory (RSPT) [47], which via the Bloch equation can be used to derive various compu- 
tational schemes, such as the linked-diagram expansion (LDE) [48, 49, 50]. A particularly powerful 
technique is the Coupled-Cluster Approach (CCA) [51, 52, 53], which is widely used in quantum chem- 
istry [54, 55]. This technique is non-perturbative but closely connected to MBPT, and we shall include 
it in the MBPT category here. The MBPT techniques are primarily developed for the weak-coupling 
case, but might in the non-perturbative (CCA) form be used also in strong coupling. 

The MBPT procedures, based initially upon RSPT, have the great advantage compared to tech- 
niques based upon BWPT that they are size -extensive at each order, which implies that the energy scales 
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Fig. 1. Schematic illustration of the connection between various many-body techniques. The upper right part 
represents the Green's-function (GF) approach, which is used to derive the Bethe-Salpeter (BS) equation, nor- 
mally analyzed in terms of the Brillouin Wigner perturbation theory (BWPT). The lower-left part illustrates the 
many-body perturbation theory (MBPT), originating from Rayleigh-Schrodinger perturbation theory (RSPT). The 
combination of quantum-electrodynamics (QED) with MBPT is represented by the covariant-evolution-operator 
(CovEvOp) method, and the link to the BS equation and the corresponding Bloch equation — the main subject of 
the present paper — is illustrated by the arrows. 



linearly with the size of the system — a property of vital importance for molecular problems [56, 57]. 
Both procedures can also be combined with the extended-model-space technique, which is particularly 
effective in dealing with problems of quasi-degeneracy [58, 59, 47, 60]. 

For QED problems the S-matrix technique has been the standard procedure since the days of Feyn- 
man and Dyson. (For a review of the application to bound-state problems, see ref. [61].) Being based 
upon scattering theory, this technique has the disadvantage that its structure is quite different from that 
of MBPT, which makes it hard to combine the procedures (see, e.g. ref. [62]). The standard proce- 
dure for such a combination has been to perform a separate (relativistic) many-body calculation and 
adding first-order energy corrections from QED analytically [63]. This procedure gives in many cases 
satisfactory results but is hard to improve in any systematic way. In particular, it gives no additional 
information about the wave function. 

Another disadvantage with the S'-matrix formalism is that the energy is conserved between the 
initial and the final states. This implies that it cannot be combined with the extended-model-space 
technique, successfully applied in MBPT. This techniques requires generally elements of the effective 
interaction that are nondiagonal in energy. This problem has recently been remedied by means of a new 
technique, known as the Covariant-Evolution-Operator method (CovEvOp), which is a modification of 
the standard evolution-operator technique of time-dependent perturbation theory [31] in order to make 
it applicable to relativistic problems (for a review, see ref. [64]). This technique has a structure that is 
very akin to that of MBPT, and it deals with the key ingredients of MBPT — the wave operator and the 
effective interaction. At the same time the method is closely related to the S'-matrix formalism and the 
Green's-function procedure. The technique can therefore be regarded as a merger of MBPT/CCA and 
QED [65], and it has recently been successfully applied to the quasi-degenerate fine-structure states of 
heliumlike systems [66]. 

The quasi-degenerate problem can also be handled with the two-times Green 's-function approach, 
developed by Shabaev and coworkers (for a review, see ref. [67]). This technique, however, has no 
direct link to MBPT and will therefore not be discussed further here. 

The procedure with the Covariant-Evolution-Operator method is now being further developed at 
our laboratory in order to combine QED and MBPT in a more complete fashion. This will be based 
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on the non-perturbative coupled-cluster approach (CCA) of electron correlation or the so-called Dirac- 
Coulomb approximation, corresponding to the "ladder approximation" of the Bethe-Salpeter equation. 
This is combined with a perturbative expansion of the remaining (mainly QED) effects, which in prin- 
ciple leads to the full BS equation. This is along the lines early drawn by Sucher [38] and followed by 
many later works [40, 15, 43, 36, 68]. Our approach differ from all the earlier ones in the sense that all 
effects are evaluated numerically rather than analytically. 

Our approach implies that the QED effects are evaluated with highly correlated (relativistic) wave 
functions, and for two-electron systems the results will then, in principle, be comparable to those of 
Drake's unified method, with the difference that the relativistic effects are included in a complete way 
and that the QED effects are evaluated numerically. 

In the diagram in Fig. 1 we have tried to represent the relations between the many-body approaches 
described here in a simple and illustrative way. The many-body procedures based upon Rayleigh- 
Schrodinger perturbation theory are indicated in the lower- left part and the Green's-function and Bethe- 
Salpeter procedures, more associated to Brillouin-Wigner perturbation theory, in the upper-right part. 
The present paper deals particularly with the connection between the two approaches, represented by 
the arrows in the diagram. 

In addition to deeper insight into the different procedures, the present treatment will make it pos- 
sible to analyze a problem based on the BS equation in terms of RS-MBPT — not only in terms of 
BWPT, as has previously been the case [40, 44]. The Bloch equation in commutator form, compatible 
with the BS equation, which is derived, has the same relation to the BS equation as has the standard 
Bloch equation to the ordinary Schrodinger equation, and it could possibly be used to eliminate the 
quasi-degeneracy problem that might appear when the BS equation is treated for a single state at a 
time. 

Since the equivalence of the MBPT-QED-CovEvOp procedure with the BS equation has now 
been established for two-electron systems, this new link will probably make it easier to apply the BS 
procedure — or its equivalence — also to systems with more electrons. Alternatively, this can be used 
to analyze a many-body-QED calculation to find out what is missing in order to represent a complete 
Bethe-Salpeter treatment. Our main emphasize here is applications to atoms and other weak-interacting 
systems. Since the procedure we have developed, however, is based upon a combination of perturbative 
and non-perturbative approaches, the results obtained might be useful also outside this regime. 

The paper will be organized in the following way. Below we shall first conjecture the Bethe-Salpeter 
equation in a simple-minded way as an introduction. In section 2 we shall summarize the necessary 
ingredients of time-independent and time-dependent perturbation theory and in the following section 
briefly review the original derivations of the Bethe-Salpeter equation by Bethe and Salpeter and by 
Gell-Mann and Low, based on Green's functions. The main part of the paper will be devoted to a rigor- 
ous derivation of the Bethe-Salpeter equation, starting from the covariant-evolution-operator method. 
The basics of the method are summarized in section 4, and the method will then be used to derive 
the Bethe-Salpeter equation. A corresponding Bloch equation will also be derived, which will make 
it possible to treat the BS equation perturbatively also for a quasi-degenerate (extended) model space. 
Technical details of the treatment are given in a number of appendices. Radiative effects (self ener- 
gies and vacuum polarization) are not considered here but can be included by modifying the electron 
propagator and photon interactions, as discussed, for instance, by Douglas and Kroll [40]. 

1.2. Bethe-Salpeter equation 

An equation of BS type can be conjectured in a very simple way by considering the time-independent 
nonrelativistic Schrodinger equation 

H^ = E^ (1) 
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with H = Ho+Vi, where Hq = hi+h 2 is the zeroth-orderHamiltonian (sum of single-electron Hamil- 
tonians) and V\ — e 2 /ri 2 is the electron-electron interaction (in relativistic units l ). The Schrodinger 
equation can then be expressed 

(E - H )V = Vi* (2) 



with the solution 
where 



# = r(£)Vitf , 0) 

^ ^ E — Hq E — £ r — £ s 

is the "resolvent" operator [47, Ch. 9] and \rs) is the Dirac notation of the straight (not antisym- 
metrized) product of two single-electron functions, satisfying the Dirac equation 

h\i) = £i\i) (5) 

We apply the summation convention, implying summation over repeated indices appearing on one side 
of the equation. Unless specifies otherwise, the summation is performed over positive- (particle) as 
well as negative-energy (hole) states. 

In the relativistic formalism one should, following Sucher [37, 38], replace V\ by A ++ e 2 /ri2A ++ , 
where A ++ is the projection operator for particle (positive-energy) states. This leads to the Coulomb- 
ladder approximation, mentioned above, i.e., a series of Coulomb interactions separated by particle 
states. In QED V\ can in the first approximation be replaced by the energy-dependent interaction with 
a fully covariant photon V\ (E), i.e., Coulomb and transverse photon, the latter representing (retarded) 
Breit interaction. In the next step V\(E) can be replaced by Vi(E) + V 2 (E), where V 2 (E) represents 
the non-separable (irreducible) interaction of two photons, i.e., the interaction of two covariant photons 
that in the QED description cannot be represented by repeated single-photon interactions (see Fig. 6 
below). Continuing this process, summing all non-separable interactions with one, two, ... photons 

V{E) = V 1 (E) + V 2 (E) + --- (6) 

leads to 

* = T(E) V(£7)tf (7) 

or 

(E - JT ) * = V(£7)tf (8) 

This is equivalent to the Schrodinger-like form of the Bethe-Salpeter equation derived by Sucher [38, 
Eq. 1.47] and also used by Douglas and Kroll [40, 3.26] and by Zhang [43, Eq. 15]. 

The BS equation (8) can be expanded in terms of a Brillouin-Wigner perturbation series [47, Ch. 

9] 

* = *o+ (r Q {E)V{E) + Tq{E)V(E)T q {E)V(E) + ■ ■ -^o (9) 
where ^ is the unperturbed wave function and 

Tq{e) = e~W q (10) 

is the "reduced" resolvent (4) with the unperturbed state removed. For this sequence to converge prop- 
erly, it is required that there be no eigenstate of H close in energy to that of * and of the same 
symmetry. A rigorous derivation of the equation will be given in the following sections. 

1 In this article relativistic units are used, i.e., m = c = h = €q = 1, e 2 = 4ira, where a is the fine-structure constant. 
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2. Conventional many-body perturbation theory 
2.1. Time-independent perturbation theory 

In time-independent many-body perturbation theory (MBPT) (see, e.g., ref. [47]) the aim is to solve 
the Schrodinger equation by successive approximations for a number of "target" states 

H9 a (x) =E a <f a (x): (a = 1,2, •••<*) (11) 

(x stands here for all space coordinates). The time-independent Hamiltonian is partitioned into a 
zeroth-order Hamiltonian and a perturbation 

H = H + H' (12) 

For each target state ^> a (x) there exists a model state or zeroth-order wave function (ZOWF) ^q(x) 
that is confined to a subspace, the model space (P), spanned by eigenfunctions of H . The model 
space can be degenerate or non-degenerate (quasi-degenerate). In the latter case the model states are 
not necessarily eigenstates of Hq. It is always assumed that all degenerate states of Hq are either 
entirely inside or entirely outside the model space. 

A wave operator £1 can be defined so that it transfers all model states to the corresponding target 
states i 

(13) 



^ a (x) = fifS(iE) ; (a = 1,2, •••<*) 



In the following we shall use the intermediate normalization (IN), implying that 

(*?(x)|* Q (x)> = 1 (14) 
The model states are the projections of the target states on the model space 

tffffc) = P^ a {x) (15) 

which implies 

PnP = P (16) 
The exact energies as well as the model states are obtained by solving the secular equation 

H eS *%(x)=E a *Z(x), (17) 

within the model space. Here, H e g is the effective Hamiltonian, in IN given by 

H cS = PHQP (18) 

The wave operator satisfies the generalized Block equation [58, 47] 



n, H ] p = (H'n - n H' eS ) p 



(19a) 



where H' cS is the effective interaction (in IN) 

H' cS = H cB - PH P = PH'ilP (19b) 

For a degenerate model space with the energy Eq the equation goes over into the original Bloch equa- 
tion [69, 70] 

(E a -H a )np= (H'n-nH' cS )p (i9c) 

The Bloch equation contains generally the information of a system of Schrodinger equations (11), 
corresponding to a number of target states. The equation can conveniently be used as the starting point 
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for generating various perturbative and non-perturbative schemes [58, 47]. It leads directly to a gen- 
eralized form of the Rayleigh-Schrodinger perturbation expansion, and it can be used to generate the 
linked-diagram expansion (LDE) as well as the non-perturbative coupled-cluster approach (CCA). The 
commutator form of the Bloch equation (19a) makes it possible to work with a non-degenerate or "ex- 
tended" model space", which is of particular importance for quasi-degenerate problems, as mentioned 
above. 



2.2. Time-dependent perturbation theory 

In time-dependent perturbation theory we start from the time-dependent Schrodinger equation 

i^- t x(t,x) = H(t) X (t,x) (20) 

As before, x stands for all space coordinates, while t is a single time variable. Even if the Hamiltonian 
may be formally time-dependent, we are interested in states that are stationary, which implies that the 
wave function has the form 

X (t,x) = <$(x)e- iEt (21) 

where E is the energy of the system and ty(x) is the time-independent wave function. The latter is then 
a solution the time-independent Schrodinger equation (2) 

HV(x) = EV(x) (22) 

In the interaction picture (IP) [31] with the partitioning (12) the wave function is related to that of 
the Schrodinger picture by 

Xi(t,x) = e iH ° t X (t,x) (23) 
and the time-dependent Schrodinger equation becomes 

i^- t Xi(t,x)=H{(t) X i(t,x) (24) 
The time-evolution operator, defined by 

Xi(t,x) = U I (t,t )xi(to,x) (25) 

then satisfies the equation 



i^U I {t,t ) = H{(t)U I (t,to) (26) 



with the solution [31, Eq. 6.23] 



Ui(t,to) = 1 + ^ i" / d 4 x„---/ dVT D [W{(a;„)W{(a;„_i)"-W{(a;i)] (27) 

n=l U - Jt ° Jt ° 

Here, x = (t, x), Td is the Dyson time-ordering operator, and TL[ (x) is the perturbation density defined 
by 

H[{t) = j ^xH'^t^x) (28) 
In applying this formalism to perturbation theory, an adiabatic damping is added [31] 

H{(t) -> H{ 7 = H{ e -Tl*l ; Ui(t, to) -» U^(t, to) (29) 
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where 7 is a small, positive number. This implies that as t — > —00 the eigenfunctions of H tend to 
eigenfunctions of H n . 

In QED the perturbation density due to the interaction between the electrons and the photon field 
is given by [71] 

n{(x) = -e^ I {x)a^A^x)i> I (x) (30) 

where e is the absolute value of the electronic charge, ^\{x), tpi(x) are the electron-field operators in 
the interaction picture, the photon-field operator and a M are related to the standard Dirac matrices 
by = (1, a). 

3. Green's function approach 

In this section we shall essentially reproduce the derivation of the BS equation by Gell-Mann and Low, 
starting from Green's functions. We consider a two-particle system for which the Green's function is 
defined [31, pp. 64 and 116] 

G(x 1 ,x 2 ;x 10 , X20) = -- — ! . . . ! — - (31) 

(Uh|Uh/ 

Here, |0h) represents the vacuum state and V> H (x), ^n{x) the electron-field operators, all in the 
Heisenberg representation. The latter are related to those in the interaction picture by 

^(t,x) = U(0,t)Mt,x)U(t,0) (32) 

where U is the evolution operator (27). Transforming the Green's function to the interaction picture 
then yields [31, Eq. 8.9], [2, Eq. 16], [64, Eq. 259] 

nl , , , (OilTp^C^^C^)^!^, -00)^ (2:20)^ (a;io)]|0i) 

G(x 1 ,x 2 ;x 10 ,X2o) = ! 7777777 77777 — (33) 

(0i|Z7i(oo, -oo)|0i) 

Obviously, only fully contracted terms contribute to the vacuum expectation value. By applying Wick's 
theorem [31, p. 83] [47, Sect. 11.5], this can be represented in terms of Feynman diagrams. The de- 
nominator has the effect of eliminating the singularities of the numerator, in the Feynman picture 
represented by unlinked or disconnected diagrams, leading to [31, Eq. 9.5] 

G{x' 1 ,x' 2 ;x w ,X2o) = -(0i|Td[^i(x'J^i(x' 2 )[/i(oo, -oo)^ I t (a; 2 o)^ I t (a;io)|Oi) conn (34) 

In contrast to the evolution operator (25), the Green's function is relativistically covariant in the 
sense that the integrations are performed over all space and time and the electron-field operators can 
represent particle (positive-energy) as well as hole (negative-energy) states. This also implies that, in 
the energy representation (fourier transform), the energy is conserved at all diagram vertices. 

The Green's function can be expressed 

G{x' ll x' 2 ]x Wl x 2 o) = G (x' 1 ,x 2 ;x w ,X2o) + 

d 4 xi d 4 x 2 d 4 x 3 d 4 x 4 G (x[, x' 2 ; xi, x 2 ) K,(xi, x 2 ; x 3 , x 4 ) G {x 3 , x 4 ] x w , ^20) (35) 

where /C represents the interaction kernel of all connected diagrams and Go is the zeroth-order Green's 
function 

G (x' 1 ,x' 2 ;xio,X2o) = -(Oi|T D [^i(x' 1 )V>i(x 2 )^ I t (x 2 o)V' I t (xio)]|Oi) 

= S , F (x' 1 ,xi )S'F(a; 2 ,a;2o) (36) 
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with Sf being the Feynman electron propagator or zeroth-order single-electron Green's function, de- 
fined by 

iS F (x',x ) = (Oi|T D [^i(a/)$(zo)]|Oi) (37) 

assuming the vacuum state be normalized. This is illustrated in Fig. 2. In operator form the Green's 
function can be expressed 

G = G + G /CG (38) 
In some cases the kernel of the Green's function can be separated into two kernels 

K = /C 2 G /Ci (39) 

with no photon-field contractions between them. The kernel is then said to be separable. If a kernel 
cannot be separated further in this way, it is said to be non-separable 2 . The complete kernel can then 
be expressed 

/C = k + kGqk + kGqkGqk + ■ ■ ■ (40) 
where k represents all non-separable kernels. This leads to the Dyson equation for the Green's function 

G = G + G q kG (41) 

illustrated in Fig. 3. 




Go 



X20 



Fig. 2. Graphical representation of the two-particle Green's function (38). K, represents all interactions between 
the electrons. 




Fig. 3. Graphical representation of the Dyson equation (41) for the two-particle Green's function, n represents the 
non-separable interactions between the electons. 



Bethe and Salpeter as well as Gell-Mann and Low argue that a related equation can be set up for the 
two-electron bound-state wave function. In that case the first (inhomogeneous) term on the rhs does not 



2 What we here refer to as "separable" and non-separable" are often referred to as "reducible" and "irreducible". Since the 
latter terms have recently been used also with a different interpretation, we avoid them here. 
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contribute, since that is in their formulation composed of free-electron propagators, and the bound-state 
wave function does not have any such components. This leads to the homogeneous equation 

* fail #2) = j j j j d 4 a;id 4 a;2 d 4 x 3 d 4 x 4 Gq{x' 1 , x' 2 ; x\, £2) k(xi, X2; X3, x&) ^(x3, X4) 

(42) 

or in short-hand notations 

* = G k * (43) 

This is the original form of the Bethe-Salpeter equation [1, Eq. 11a], [2, Eq. 37]. It should be noted that 
this wave function contains individual times for the two particles. This reflects one of the problems 
referred to in the Introduction. The relative time between the particles does not correspond to any 
physical quantity and leads to spurious solutions. There are several ways of eliminating the extra time 
dependence in a covariant way. Sucher [38], following Salpeter [18], integrates the fourier transform 
over the relative energy, which leads to a Schrodinger-like form with a single time/energy dependence 
of the type (8) given above. This reduction can be done without loosing any physical content of the 
original equation [15, 16, 32, 17]. In the following sections we shall derive an equivalent equation in a 
different way. 

Our notations here differ from those used by Bethe-Salpeter and Gell-Mann-Low. The Green's 
function (31) is in their works denoted by if (12, 34) and referred to as the "amplitude function for 
the propagation of the particles" by Bethe-Salpeter [1] and as the "two-body kernel" by Gell-Mann- 
Low [2, Eq. 11]. Our "non-separable kernel" k is by BS denoted by G and referred to as "irreducible 
graphs" and by GML denote by G and referred to as the "interaction function". 

4. Covariant evolution operator approach 
4.1. Definitions 

In the following sections we shall derive the Bethe-Salpeter equation, starting from the covariant form 
of the evolution operator [64]. This will demonstrate the relation between the BS equation and standard 
many-body perturbation theory (MBPT) in a clear way. In the present section we shall first review the 
basics of the evolution-operator method and in the next section use that method for deriving the BS 
equation. This will directly lead to the Schrodinger-like form (8). 

According to the Gell-Mann-Low theorem [2, 31, p. 61] the time-independent wave function (21) 
can in the case of a single target function be expressed in intermediate normalization (IN) (14) as 

*t->-**->-ft <t.p;(o,-Ui».> <44) 

where J7 7 is the evolution operator (29) and "Po ( x ) is the time-independent zeroth-order wave function 
(15). (From now on we work in the interaction picture and leave out the subscript "I".) is an 

eigenfunction of the Hamiltonian H n + H' 

(H + H')^(x) = E^(x) (45) 

where H 1 is in our case the electron-field interaction (30). Since this perturbation represents an uncon- 
tracted photon, the wave function will generally lie in an extended Fock space, where the number 
of photons is not conserved. 

The GML formula can be generalized to a general multi-dimensional model space [64, Eq. 110] 

v ' 7 ^ ($" ^ (0,-00) $° v 7 
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where the function <f>" is defined 



<fr a (x) = lim lim x a (t, x ) 

7^0 t^ — oo 



(47) 



This function is generally distinct from the zeroth-order wave function (15) in intermediate normaliza- 
tion. Since the function (47) generally does not satisfy IN, a normalization constant TV" is inserted. 



t = t' 



^Noncov(^^o) — 1 + 



t = t 



K 



Fig. 4. Graphical representation of the non-covariant evolution operator (48). The time evolution occurs only in 
the positive direction. 



t=t 




K 

4> 



t=t' 
t=i 



ft 



Fig. 5. Graphical representation of the covariant evolution (49) (left). Here, time evolution can occur in the pos- 
itive as well as the negative direction. The right part of the figure depicts the relation to the two-times Green's 
function (50). 



For a two-electron system the non-covariant evolution operator (25) can in analogy with the 
Green's function (38) be expressed 

t/Noncov(t',M = l + ft + {x[)ft + (x' 2 )lCi> + (x2o)i>+(x W ) (48) 

where again K, represents the kernel of all fully contracted (separable and non-separable) interactions 
and ft,, ip + the positive-energy part of the electron-field operators. This is illustrated in Fig. 4. In 
contrast to the Green's function above, the evolution operator (27) has a single initial time t = to and a 
single final time t = t'. The time integration is performed from t = to to t = t' — only in the positive 
direction — which implies that the operator is not relativistically covariant. 

A fully covariant form of the evolution operator that is applicable to relativistic problems can be 
obtained by inserting electron propagators in the non-covariant expression, as indicated in Fig. 5 (left), 
corresponding to the expression [66, 64, Sect. 5] 

EWi',*o) = 1 + jjjj d 3 xi d 3 x' 2 d 3 x w d 3 x 20 ft 0^)^04) G /CG feoMzio) ( 4 9) 

leaving out the integrations over the coordinates of /C (see Eq. 35). It then follows from the relation 
(38) that the covariant evolution operator is related to the two-times Green's function (where all initial 
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and all final times are equal) by 

Ucov(t',t ) = JJJJ d 3 x' 1 d 3 x' 2 d 3 x w d 3 x 20 ft (x'jft (x' 2 ) G(x[, x 2 ; x w , x 20 ) ip(x 2 o)ip(x w ) 

(50) 

as illustrated in Fig. 5 (right). 

From the relations [64, Eq. 193 (note misprints)] 

J d 3 x iS F (x,x a )i}(x ) = Q(t - t ) $+{x) - 6(i - t) ft(x) 

J d 3 xft(x)iS F (x,x a ) = e(t-t )ft + (xo)-e(t a -t)ft_(x ) (51) 

it follows directly that the form (49) is equivalent to the non-covariant form (48), when only particle 
states are involved. That the former in addition is relativistically covariant follows from the fact that the 
electron-field operators can represent particle as well as hole states and the internal time integrations 
are performed over all times — in the positive as well as the negative direction. From now on we shall 
work only with the covariant form of the evolution operator and leave out the subscript " Gov"- 

In using the evolution operator in perturbation theory, we assume that we operate to the far right 
on positive-energy states in the model space. Then, as shown in Appendix .2, we can eliminate the 
rightmost zeroth-order Green's function and set the initial time to t = — oo. We shall also assume that 
the limit of the adiabatic damping 7 — > is taken. 

The Covariant evolution operator is closely related to the Green's function — the main difference 
being that the Green's function is a function, while the evolution operator is an operator. The poles 
of the Green's function (in the energy representation) correspond to the energies of the system, while 
it gives no direct information about the wave function. The covariant evolution operator, on the other 
hand, contains information about the energy as well as the wave function. 



4.2. Model-space contributions 

Even after eliminating unlinked or disconnected contributions in Eq. (34), the evolution operator may 
contain (quasi)singularities, namely when the intermediate state of a separable kernel lies in the model 
space and is degenerate or nearly degenerate (quasi-degenerate) with the initial state. As mentioned, 
a kernel is said to be separable, if it can be separated into two kernels with no photon contractions 
between them. Singularities appear only for separable interactions. In the covariant-evolution-operator 
approach these singularities are eliminated by introducing a reduced evolution operator U(t, —00) [66, 
64, Eq. 116], defined by 

U{t, -oo)P = P + U(t, -oo)P ■ PU(0, -oo)P (52) 

Here, the last term is a product of two operators that evolve independently from an initial state in the 
model space (t — —00), which is indicated by the "dot". Note also that the last factor has the final time 
t = and hence is time independent. This situation should be distinguished from the case where two 
operators are "coupled" and operate "in succession" 

U{t,t ) = U(t,t")U{t",t ) (53) 

This distinction will be important for the following treatment. 

Normally, we shall assume that the initial time in the evolution operator is to = —00, and in cases 
where there is no risk for ambiguity we shall leave that out from the operator, so that 

U(t) = U(t, -00) 
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The definition (52) will then be written 



U(t)P = P + U(t)P ■ PU(0)P 



We also introduce the notation U'(t) = U(t) — 1, which yields in place of the definition (54) 

U(t)P = U'(t)P - U(t)P ■ PU'{0)P 
Here, the last term is the counterterm 



C(t)P = -U(t)P ■ PU'(0)P 



which removes the (quasi)singularities. This can also be expressed 

C{t)P = -U{t)P • PU(0)P - U{t)P • PU(0)P • PU{0)P 



(54) 



(55) 



(56) 



(57) 



After removing a singularity, there is normally a non-vanishing remainder, referred to as the model- 
space contribution (MSC), defined as 



MP = U(0)P - U(0)P 



(58) 



and further discussed in the Appendices. The new operator U ( "U-bar") is defined as the evolution op- 
erator with all model-space states removed. (The MSC is analogous to the reference-state contribution, 
appearing in the S'-matrix formalism, where the effect normally appears only when the intermediate 
states is equal to the reference or initial state. In our formalism with an extended model space the effect 
can appear also for other model-space states, and we prefer the more general term.) It should be noted 
that the counterterms also remove quasi-singularities, due to quasi-degenerate states that are included 
in the model space. This can be of vital importance for the convergence of the procedure. 

As discussed in Appendix .3, the model-space contributions are of two kinds. The first kind appears 
for all interactions, even if they are time or energy independent, while the second kind appears only 
for time- or energy-dependent interactions. The first kind appears also in standard time-independent 
perturbation theory and corresponds to so-called folded diagrams of MBPT [47, 64, Fig. 5]. 

4.3. The wave operator and effective interaction 

As mentioned previously, the evolution operator (27) with the perturbation (30) can contain uncon- 
nected photon operators, which implies that it operates in a general Fock space, where the number of 
virtual photons is not conserved. We then separate the covariant evolution operator (49) into 



U(t) = PU(t) + QU(t) 



(59) 



where Q = 1 — P is operating in the general Fock space, while P is the projection operator for the 
model space, confined to the restricted Hilbert space with no uncontracted photon. This leads with the 
definition (54) of the reduced evolution operator for t = to the factorization theorem [64, Eq. 121] 



U(0)P = [1 + QU(0)] P • PU(0)P 



where the first factor on rhs is regular. Inserted in the GML formula (46), this yields 

* a (a0= [l + QU(0)] ^ 



(60) 



(61) 
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where ^>q is the zeroth-order wave function (ZOWF) (15) in intermediate normalization 
The square bracket above is the wave operator 

n = i + qu(o) 

* a (x) = n^(x) (63) 

The result here is a direct consequence of the generalized Gell-Mann-Low theorem and the definition 
of the reduced evolution operator. 

As mentioned, with the perturbation (30) the wave function H> a lies generally in a Fock space 
where the number of (virtual) photons is not conserved. But we are interested here in the case where 
all photon operators are fully contracted, and for that purpose we project the equation on the restricted 
Hilbert space without uncontracted photon operators 

TV a {x) =V[1 + QU(0)] %{x) (64) 



or 



tf a (x) = [l + QI7(0)]*S(aO (65) 

where = V'9 a (x) is the projected wave function on the restricted Hilbert space and Q = VQ 

is the conventional projection operator for the complementary space (outside the model space). The 
wave operator in this space is 

Q = Vft = 1 + QU(0) (66) 



In IN (14) the wave operators satisfy in both spaces the relation (16) 

PftP = PflP = P. (67) 
The effective interaction (19b) is in this formalism given by [64, Eq. 130] 



K« = p 



P 

t=0 



(68) 



5. Connection to the Bethe-Salpeter equation 
5.1. Expansion of the wave operator 

We know from the generalized Gell-Mann-Low relation (46) that the wave function *& a (x) in the 
extended Fock space satisfies a Schrodinger-like equation (45) with the Hamiltonian H = Ho + H' , 
where H' is the perturbation (30). We now want to find the corresponding equation for the wave 
function v?" = P<P a (x) in the restricted space with no uncontracted photons, and we shall see in this 
section that this leads to the Bethe-Salpeter equation. 

We shall start with the exchange of a sequence of separable covariant single photons between the 
electrons, which can then be generalized to other interactions, leading to the full equation. This will 
first be done for a degenerate model space and then extended to the general case. 

As shown in Appendix .2 (Eq. 164), the contribution to the wave operator due the exchange of one 
single photon is 

fiWp - Qt/ (1) (0, £)P = T Q {£)V(£)P (69) 

and the corresponding contribution to the effective Hamiltonian H^g (£) = PV(£)P. Here, Tq(£ ) is 
the "reduced" resolvent (10) and V{£) is the effective single-photon potential (162), assuming that we 
operate to the right on a fourier transform (157) with the energy parameter £ . 
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Similarly, it is demonstrated in Appendix .3 (Eq. 180) that the contribution to the evolution operator 
from two separable single-photon interactions is for a degenerate model space given by (leaving out 
the arguments) 

n p)p = qu^p = T Q vnWp + ^ * h£ (70) 

where the last term represents the model-space contribution (MSC) (58) 



(The asterisk is introduced here only to indicate that there is a cancelled singularity at that position, 
which is of importance for the further treatment, as discussed in the Appendices.) The contribution to 
the effective Hamiltonian (182) due to two-photon exchange is 

H$ = PVn^P + d J§- * H$ = H% + ^ * H$ (71) 

The last term is the MSC to the effective interaction, and if the model space is degenerate with the 
energy Eq that term becomes 



9H * * P dV{£) 



PV(E )P (72) 

En 



This corresponds to the "reference-state contribution", discussed in connection with the S'-matrix treat- 
ment of two-photon exchange [72, 73]. 

The treatment above will now be generalized to all orders as a first step towards deriving the full 
BS equation. We start with the covariant evolution operator (49) U(t) = U(t, — oo) and the reduced 
evolution operator (55) 

U(t)P = U'(t)P - U(t) • PU'(0)P (73) 

where U' = U — 1. Note that only the first factor in the product is time dependent (see. Eq. 172). Note 
also the appearance of the "dots" in this expression. The significance of the dot is discussed in relation 
to the definition (52). 

In the following we shall leave out the prime on U' and also the time arguments, if there is no risk 
of ambiguity. We then express the counterterm (56) as 



CP = —U • PUP 



(74) 



and the evolution operator is given by 

UP = TVP + TVTVP + YVTVYVP + ■■■ (75) 

where T = T(£) is the resolvent (4). The "U-bar" operator (58), with all intermediate model-space 
states removed, is 

UP = TVP + TqVTqVP + TqVTqVTqVP +■■■ (76) 
We introduce a special symbol for the time derivative at time t = 



: .OA 

A = l-r- 

dt 



t=o 



(77) 



Since the evolution operator (159) has the time dependence 

U(t,£) = e-^ £ - H ^U(0,£) 
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it follows that the time derivation eliminates the denominator of the first (leftmost) resolvent, so that 
PUP = P(V + VTV + VTVTV + ■ ■ -)P = PVP + PVUP 

PUP = P(V + VT Q V + VT Q VT Q V + ...)P = PVP + PVUP (78) 
The effective interaction H' cS (68) is with this notation given by 



H' eS = PUP 



(79) 



We also introduce the corresponding "H-bar" operator with no intermediate model-space states 

H' eS = PUP (80) 
We recall the definition (58) of the model-space contribution (MSC) 

UP = UP + MP (81) 

and can easily derive the identities 

UP = UP + U PUP = UP + U PUP - M PUP (82) 
UP = UP -U PUP + U PUP PUP (83) 

Then the reduced evolution operator (73) becomes 

UP = UP + CP = UP + U PUP - U ■ PUP - M PUP (84) 
which using the definition (81) leads to the series 

UP = UP+ (U PUP - U • PUP) (l - PUP + PUP PUP +■■■) (85) 
With the identity (83) this becomes 



UP = UP+ (U PUP - U • PUP) 



(86) 



which is an exact expression also for a quasi-degenerate model space. It can be expanded as 

UP = UP+ (U PUP - U • PUP) + (U PUP - U • PUP) (PUP - -PUP) + ■■■ (87) 
As discussed in Appendix .3, the result (86) can be expressed 

UP=UP+ — * PUP = UP+ — *H' cS (88) 
ot ot 

where S£ is the change in the model-space energy, represented by the "dot", SU is the corresponding 
change in U, and H' cS is the "H-bar" operator (80). In the case of complete degeneracy this becomes 



at 



* H' e s (89) 

£=E 



Introducing the" Omega-bar" operator ft (with no intermediate model-space states) in analogy with 
the wave operator (66) 

nP = P + QUP = P + T Q VP + T Q VT Q VP+--- (90) 
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we can express the relations above as 



np = np 



S£ HeS 



on 

~d£ 



£=E a 



*H> S 



(91) 



The second term is here consequently an exact expression for the entire model-space contribution to 
the wave operator. This is in agreement with the three-photon result (191). 

By taking the time derivative of the relation (86), using the relations above, we obtain similarly 



H' eS = H' ef 



SHL 



* H' cS => H' e s 



d£ 



*H' 



-E 



(92) 



The second term represents here the model-space contribution to the effective interaction. This result 
agrees also with the third-order result (193). 

From the results above we conjecture that the wave operator can at complete degeneracy alterna- 
tively be expressed 



np 



3! d 3 £ 



^ n\ d£ n 1 eff; 



n=l 



(93) 



with all derivatives taken at £ = E , and we shall now prove this relation by showing that it is compati- 
ble with the results (91) and (92), which we have rigorously derived. This equation contains eliminated 
singularities, indicated by the asterisks. As discussed in the Appendices, the derivative of such an ex- 
pression has to be taken before the singularity is eliminated. Using the rules developed, particularly in 
Appendix .5, we find for instance 



d£ V S£ cff 



2 8 2 £ cff 



dn dHL 



d /is 2 n 

d£\2~5 2 ~£ 



cirj 



d£ d£ 



i d 3 n , , L „ ... 



dH. 



off 



d£ 



* H 'cff 



(94) 
(95) 



Note that in the second example the two H' cS operators have in the quasi-degenerate case different 
energy parameters, and therefore only one of them is affected by the derivation. 

Generalizing these rules, we can evaluate the derivative of the wave operator (93) 



d£ 



+ 



dfl 
~d£ 
dSl 
~d£ 



2&£ cfr 



i d 3 n 
3! We 



or 



d£ 

dfl 
~d£ 



= E 



n=l 



id 2 n 

2d 2 ! 

1 d n Cl 
n\ ~d£™ 



dHL 



off 



d£ 



* (H' cS ) 
*H' cS + 



1 d 3 ti dHL 



cir 



3! d 3 £ d£ 
dH' 



* (H' cS f + 



' eS > ~d£~ ^ cS ' 



(96) 



(97) 



We now insert this expression into the equation (91), which yields 



d£ tL 



off 



i d n n 



71=1 



[ H cs) 



H'cS 



d£ ™ } 



1 Hcs 



or, using the relation (92), 



n—1 



(98) 



(99) 



This is identical to the conjectured relation (91) and therefore completes the proof. The sum represents 
by definition the model-space contribution (MSC). 
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Fig. 6. Examples of non-separable two-photon interactions. 



5.2. Derivation of the Bethe-Salpeter equation. Degenerate model space. 

The previous treatment has been based upon the Hamiltonian H = Ho + V(E), where V(E) is the 
potential due to the exchange of a single covariant photon. But the process can be repeated in exactly 
the same way, if we include all non-separable multi-photon interactions. A non-separable interaction 
is defined as an interaction that cannot be represented by two or more simpler interactions in the way 
treated here. Two photons — crossing or noncrossing — that overlap in time represent non-separable 
two-photon interactions (see Fig. 6, c.f. also Ref. [1, Fig. 1]). These can also include the radiative 
self-energy and vertex corrections. In a similar way non-separable three-, four-,... photon interactions 
can be defined. Therefore, in the following we replace the single-photon potential V by the general 
potential due to all non-separable interactions 



V(E) = V(E) + V 2 (E) + V 3 (E) 



(100) 



As discussed in the Appendices, when operating on a fourier transform of definite energy £, the 
energy parameter of ClP is equal to that energy, i.e., 



QF(£) = n(£)F(£) 
For a degenerate model space of energy E this means that 

ClP = ti(E )P = P + Tq(E )V(E q ) + 



(101) 



(102) 



The model functions are eigenfunctions of the effective Hamiltonian (17), and for a degenerate 
model space (of energy E ) they are eigenfunctions also of the effective interaction (19b), 



(103) 



where AE a = E a — E - Operating with the operator equation (99) directly on the model function 
then leads to the Taylor expansion 



n(E ) + - 



e 2 8 2 £ 



(AE C 



1 

+ 3T 



% = n(E a )% 



(104) 



This implies that the MSC term shifts the energy parameter of the resolvent as well as that of the 
potential from the unperturbed energy E to the exact energy E a . But Cl(E a ) with the energy 
parameter equal to the full energy for the state ^ a is also identical to the Brillouin-Wigner expansion 
(9), 



n(E a )^ = 



i 



_Q 

E a - H 



V(E a ) 



_Q 

E a ~H Q 



V{E a ) 



E° 



H< 



V{E a ) + 



o (105) 
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which represents the full wave function, i.e., 

n(E a )% = (106) 



This implies that the relation (99) essentially represents the link between the Rayleigh-Schrddinger 
and the Brillouin-Wigner expansions for an energy-dependent interaction and at the same time the link 
between the MBPT approaches and the Bethe-Salpeter equation (indicated by the arrow in the diagram 
of Fig. 1). 

The BW expansion (105) can be expressed 

ntfg = % + r Q (E a )v(E a )n^ (107) 

or 

(E a - H ) Q*" = QV(E a )^ a (108) 
From the relation (92) it can be shown in analogy with the relation (99) 

With the definitions (80) and (90) this leads to 

H^ = PV{E )n{E )P (110) 
and in analogy with the relation (104) to 

H' eS = PV{E a )Cl{E a )P = PV(E a )nP (111) 
This leads together with Eq. (108) to the final equation 



(E a - H )^ a = V(E a )^ a (112) 



This is the Bethe-Salpeter equation for energy-dependent interactions in the Schrodinger-like 
form (8). 

We have now confirmed that the Schrodinger equation (45), obtained directly from the general- 
ized Gell-Mann-Low relation in the extended Fock space with the perturbation (30), corresponds in 
the projected Hilbert space with no uncontracted photons to a Schrodinger-like equation with the per- 
turbation (100). Both forms represent the complete interaction between the particles and are exactly 
equivalent to the original Bethe-Salpeter equation (42). 

The main difference between the original form of the BS equation and the Schrodinger-like form 
derived here is primarily that the latter has the time dependence reduced to a single time, which makes 
the wave function in accord with standard quantum mechanics. Furthermore, the Schrodinger-like form 
contains explicitly the resolvent, while the remaining part of the Green's function (139) is merged with 
the kernel k to form the potential V. 

The Schrodinger-like equation (1 12) we have derived is equivalent to the equation derived from the 
BS equation by Sucher [38, Eq. 1.47] and rederived by Douglas and Kroll [40, Eq. 3.26] and Zhang [43, 
Eq. 15]. In these works the equation is essentially obtained by integrating over the relative energy of 
the particles, thereby transforming the equation to an "equal-times" equation. This equation is then 
analyzed in terms of the Brillouin-Wigner perturbation theory. In our presentation the corresponding 
equation is obtained by starting from MBPT in the Rayleigh-Schrodinger formulation and summing 
all relevant perturbations to all orders. The present derivation therefore can serve as a link between the 
two approaches. 

In the next section we shall extend the treatment to the quasi-degenerate case and derive the corre- 
sponding Bloch equation. 
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5.3. Derivation of the Bethe-Salpeter-Bloch equation. Quasi-degenerate model space. 

We have previously assumed that the model space is degenerate, which for a two-electron system 
implies that the effective interaction is diagonal within this space (assuming the basis functions have 
definite symmetry). Then the relation (103) simplifies the treatment, and the formulas derived in the 
previous section lead directly to the standard Bethe-Salpeter equation (112). The treatment above, 
however, is more general and can be extended to the case where the model space is non-degenerate 
(quasi-degenerate). In the present section we shall show how this can be performed. 
The following relation can easily derived by induction 

W- TQ ~d£^~ nTQ d£W (U3) 
To prove this we form the next-order derivative 

dS^+V ~ Q d£ n + Q d£( n+1 ~> + Q d£( n -V ~ Q d£™ 

(Since no singularities are involved here, ordinary rules of derivation can be used.) Inserting the ex- 
pression (1 13) in the first term, yields 

In first order we have with = 1 — TqVQ 

on djvti) == 
d£ ~ rQ ^£-~ TQn 

which completes the proof of the relation (113). 

The formula above leads together with the expansion (99) to 

Q n P = Q^p + r Q £ 1 ^ (K sT - r «E J^~iy. W^T ™" (1 14) 

n=l n=l v ' 

The first term on the rhs can also be expressed TqV&P, and the last term is simply — TqCIH^ s , which 
yields 

qnp = r Q vnp + r Q ]T ^ -^r 1 (K s ) n r Q n^ ff (lis) 

n=l 

We can consider VO as a single energy-dependent operator, and if that operates on a particular 
model state of a degenerate model space of energy Eq, the first two terms of the bracket above repre- 
sents the Taylor expansion 

V(E )n(E ) + y2--±—l (AE a ) n = V(E a )n(E a ). (116) 

Tl. UO En 

n=l 

Thus, the expansion has the effect of transforming the energy parameter of the product Vf2 from Eq to 
the full energy E a , 

V(EoME )^o -» V(E a )n(E a )V° (117) 
in analogy with the expansion (105). Using the relation (106), the equation (115) above then becomes 

Qfittff =T Q [V(E a )n-nH' cfi ]y° (118) 
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or 

(Eo - H Q )Q^ =Q[V(E a )Q - fifl^]*? (119) 

which is consistent with the Bethe-Salpeter equation (1 12). 

If the model space is non-degenerate (quasi-degenerate), then the relation (103) is no longer valid, 
and the expansion (115) can not be expressed by means of a single energy parameter as in the Taylor 
expansion (116). Instead, the potential will depend on the full matrix of the effective Hamiltonian. We 
then replace the energy parameter in (102) by the model Hamiltonian H , 

QP = P + T Q (H )V(H a ) + ■■■ = Q(H )P 

By this notation we understand — in accordance with the rule (101) — 

A(H )B<f> = A(E )B^ (120) 

when $ is an eigenfunction of Ho with the eigenvalue Eq and B is an arbitrary operator combination. 
Together with the linearity condition, 

A(H )B(a<S> + &$') = aA(E Q )B<S> + bA(E' Q )B<S>' (121) 

where $' is another eigenfunction of Ho with the eigenvalue E Q , this defines the notation fully. 

The expansion (115) can now be regarded, in analogy with the energy modification (116), as mod- 
ifying the parameter H a to the full effective Hamiltonian H cS = H + H' cS 



n— 1 

i.e., 

V(H )ti(H )P -> V(H eS )Cl(H eS )P (123) 



and Eq. (118) becomes 

QilP = r Q (H )V(H cS )n(H cB )P - T Q (H )nH' eB (124) 
The notation here is defined by the relation 

A{H eS )B% = A(E a )B^ (125) 

where is a model function (eigenfunction of H e s with the eigenvalue E a , see Eq. 17), which 
together with the linearity condition defines the operator when acting on any model space. 
Similarly, the expansion (99) yields 

QP = ti(Ho)P + Y.~^ W = ^eff)P (126) 

n=l 

and we can now express the equation (124) as 

QfiP = r Q (H )V(H eS ) flP - r Q (H )ClH' eS (127) 

or 

TgiHo^QQP = Q[V(H eS )Q- QnH' cS ]P (128) 
Operating on a model-state function ^(Eg ) of energy Eg , we have according to the definitions above 

[TqIHo)}- 1 ^) = (E£ - Ho) n$(EZ) 
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Therefore, the inverse of the resolvent can be expressed as a commutator 

[T Q {H )]- l A^[A,H ] 
where A is an arbitrary operator. This leads to the commutator relation 

[SI, H }P = QV(H cS ) np - Qft H' eS 
The relation (111) can be generalized to 

h' cS = pv(h cB ) np 

and with the IN relation (67) PflP = Fwe arrive at the BS equation in commutator form 



Q,H ]P = V(H eS )QP-QH^ 



(129) 



(130) 



(131) 



(132) 



The equation above — which we shall refer to as the Bethe-Salpeter-Bloch equation — is the main 
result of the present work. It has the same relation to the standard BS equation as has the standard 
Bloch equation (19a) to the ordinary Schrddinger equation. It can be used to generate the BS equation 
perturbatively, essentially as the ordinary Bloch eqution is used in standard MBPT. The commutator 
form makes it possible to apply the equation to an extended model space, which essentially eliminates 
the quasi-degeneracy problem that might appear in applying the standard equation directly on a single 
state. 



5.4. Expansion of the Bethe-Salpeter-Bloch equation 

I 



I 



I 




(a) 



(b) 



(c) 



(d) 



Fig. 7. Illustration of the perturbative-nonperturbative procedure for solving the Bethe-Salpeter equation, de- 
scribed in the text. 



The original Bethe-Salpeter equation contains the exact energy and is therefore normally treated 
by means a Brillouin-Wigner perturbation expansion [40, 44], which requires a self-consistent treat- 
ment. The BS equation in the Bloch-equation form can be used to generate a perturbation expansion of 
Rayleigh-Schrodinger type that does not require any self-consistence procedure. We have at our lab- 
oratory developed a procedure of solving the Bethe-Salpeter-Bloch equation that is a combination of 
perturbative and non-perturbative techniques, which we shall here briefly indicate. A detailed descrip- 
tion of the procedure together with numerical results will appear shortly [74]. 

Our procedure is based upon the iterative solution of pair equations [75, 76, 77, 78, 47]. This 
represents the "ladder" approximation of the BS equation, as indicated in Fig. 7 (a). The pair function 



NRC Canada 



Lindgren 



23 



is combined with the emission of a single (uncontracted) photon (Fig. 7 b). This represents a function 
in the extended Fock space, discussed in section 4.3. This function can be iterated further, before 
the photon is annihilated, which can yield instantaneous (Coulomb and Breit) interactions, crossing 
the photon. These iterations can be continued after the annihilation, as indicated in Fig. 7 (c). By 
annihilating the photon on the same electron line, leads to self-energy and vertex corrections (Fig. 7 
d). At present time it is possible to treat only one co variant photon in this way, but the dominating part 
of the multi-photon exchange will be included by the crossings with the instantaneous interactions (c). 
This will correspond to all effects treated by Zhang [44] in his analysis of the helium fine structure up 
to order ma 7 , except for the non-separable part of two-photon exchange (Fig. 6). With our approach 
this part has at present to be included analytically. 

6. Summary and conclusions 

Standard many-body perturbation theory (MBPT) is conveniently based upon the Bloch equation, 
which is the generating equation for Rayleigh-Schrodinger perturbation expansion. The Bloch equation 
can also be used to generate various other perturbative schemes, such as the linked-diagram expansion, 
and it also leads to non-perturbative schemes, such as the Coupled-Cluster Approach. In the commu- 
tator form (19a) the Bloch equation leads to schemes that can handle the quasi-degenerate problem in 
an efficient way by means of an "extended" model space. 

In this paper we have reviewed the connection between relativistic MBPT and quantum-electrodynamics 
(QED) for a two-electron system by means of the recently introduced covariant-evolution-operator 
method [64]. The exchange of a single covariant photon is treated to all orders, and this is shown to 
lead to an equation of the Bethe-Salpeter (BS) type. Extending the treatment to all non-separable inter- 
actions (including radiative corrections) leads to the full BS equation. This establishes a link between 
the perturbative and non-perturbative schemes, based upon Rayleigh-Schrodinger perturbation theory 
and schemes based upon the BS equation, which are normally treated by means of Brillouin-Wigner 
perturbation theory. 

In addition, a Bloch equation in commutator form that is compatible with the BS equation is de- 
rived. This equation has the same relation to the Bethe-Salpeter equation as has the standard Bloch 
equation to the ordinary Schrodinger equation and represents a series of BS equations, associated with 
a model space that need not be degenerate. This can be used to generate a perturbative expansion, 
corresponding to the BS equation for an extended model space. In principle, this will make it possible 
to treat the quasi-degeneracy problem also within the BS formalism. Such a scheme is presently being 
tested at our laboratory. 
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UJ\ = UJ 



r s i 



UJ 2 = £ — W 



ho 



ho 



Fig. 8. Graphical representation of the zeroth-order Green's function (133). 



APPENDIX 



.1. Zeroth-order Green's function 

The zeroth-order Green's function (36) in Fig. 8 is in coordinate representation 

G (xix 2 ; X10X20) = S F (xi;x w ) S F (x 2 ; x 20 ) 
where Sf is the electron propagator 



S(x,x ) 



2tt 



S(u>) e 



-iui(t— to) 



with the fourier transform 



(x\S(lj)\x ) = = (x\S(u)\x ) 



UJ — e r + lT]r 



and the corresponding operator form 



S(uj) = 



A+ 



+ 



A_ 



uj — h + 'vq uj — h — 17] 



(133) 



(134) 



(135) 



(136) 



Here, h is the single-electron Dirac Hamiltonian in the field of the nucleus and A± are projection 
operators for positive and negative-energy single-particle states. 

We consider the equal-times Green's function with h=h = t, which gives 



Go(t,Xl,X2-,Xlo,X2o) = J ^ 

duj r 1 
2^ 



{x 1 x 2 \rs)(rs\x w x 2 o) 



€ £ r £ g 



iuj — e r + ii] r e — uj — e s + vq s 
with x = (t, x), Ldi = uj and e = uji + uj 2 . The fourier transform with respect to t is then 



(137) 



Go{e,xi,x 2 :x w ,X2o] 
(x 1 x 2 1 rs)(rs \ x 10^20} f duj 



6 £ r £ s 

or in operator form 



r dw 

J 2ir 



1 



+ 



1 



uj - e r + ir] r e - uj - e s + ir] s 



duj 



Go(e)=T(e) / ^-g (e,uj)e^ tw e^'^ 



(138) 
(139) 
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where T(£) is the resolvent (4) 



£-H 

Hq = h\ + h 2 is the zeroth-order Hamiltonian (30) and 



(140) 



g (e,u) = A+ 
+ A_ 



— + 



1 



ui — hi + it] e — u — h 2 + ii] 
1 1 



.w — hi — vq e — uj — hi — \r\ 
The inverse transformation is 





G (i,a;i,a;2;ario,a;2o) = J 


de 


e- ie * G (e,SBi 


x 2 ;x w ,x 20 ) 


and specifically, 












Go(t = 0,xi,x 2 :x w ,x2o) 


-\ 


f 2~^ G °( e,Xl > 


x 2 ;x w ,x 20 ) 



(141) 



(142) 



(143) 



.2. Single-photon exchange 

(See Ref. [64, Eq. 312].) We consider now the covariant evolution operator (49) for the exchange of a 
single covariant photon, represented by the diagram in Fig. 9 (left) 

U^{t',t a ) = - X - ff d 3 x'id 3 x' 2 ^(x'i)^(x' 2 ) ff d 4 xid 4 x 2 iS F (x'i,xi)iS F (x' 2 ,x 2 )iI(x 2 ,xi) 
d 3 x w d 3 x 20 iS F (xi,x w ) iS F (x 2 ,x 20 ) ip(x 2 o)ip(x w ) (144) 



leaving out the damping factors. More compactly, we express this as 

UW(t',t ) = -^(x' 1 )^(x / 2 )G a (x' 1 ,x' 2 ;x 1 x 2 )iI(x 2 ,xi) 
x G (xi,x 2 ;x w x 20 ) ■tp(x 2 o)'4'(xio) 



(145) 



with integrations over all variables that do not appear on the left-hand side. Here, I(x 2 ,xi) represents 
the single-photon exchange 



I{x 2 ,xi)=ea'iD FtiV {x 2 -xi)ea v 2 = J ^ e^* 2 "* 1 ) il(z, x 2 , xi) 



(146) 



where D Ftlv (x 2 — xi) is the Feynman photon propagator. 

If we operate to the right on a positive-energy state, we can use the relations (51) to simplify the 
expression. Furthermore, since in that case to < ii, t 2 and since t\,t 2 run from — oo to +oo, we must 
have to = — oo, yielding 

U^{t', -oo) = ^(x / 1 )^(x' 2 )G (x' 1 ,x' 2 ;x 1 x 2 )iI(x 2 ,xi) i>(x 2 )i>( Xl ) (147) 



The electron-field operator is in the interaction picture [31] 

$(x) = i>{t,x) = Cj4>j(x) e~ 1£jt 



(148) 
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t = t' --- 



t = t --- 




i = t' -+ 

1 



r s 



cj' = £ — uj 
2 

LU 2 = £ — Wl 



Fig. 9. Graphical representation of the covariant-evolution operator for single-photon exchange in the form (144) 
(left) and in the form (147) with to — ► — oo (right). 



with the fourier transform 

$(u},x) = j &te lu}t -4>{x) =c j (p j (x)2Tr5{uj-e j ) (149) 

(as usual, summed over repeated indices) and the inverse transform 

4(x)= J ^e-^^x) (150) 

An arbitrary function of x\ = (t\,x\) can be fourier expanded as 

F{ Xl ) = J ^e^% lJ;i ) (151) 

Operating on a (time-independent) fourier component of that function with the electron-field operator 
(150), yields 

/j 
— e- i " t i>(u>,x)F(u, 1 ,x 1 )=e- i ^ t F(u> 1 ,x 1 ) (152) 

With the adiabatic damping the time-independent component corresponds to t\ = — oo, which implies 
that the field operator propagates the function F(wi, X\) from the time — oo to t. Similarly, the product 
of two electron-field operators tp(xi)tp(x 2 ) operating on a time-independent two-electron function, 
propagates the individual electrons from the time t = — oo to t = t\ and t = t 2 , respectively, without 
any electron-electron interaction. Thus, 

j>(xi)$(x 2 ) F(w u w 2 ) = e- itlWl e- it2 ^ F(u> u uj 2 ) (153) 

We now use the form (142) of the Green's function 

G (x[,x' 2 ;x 1 ,x 2 ) = J ^e- kt 'G (e) = J ^ e~^' T(e) J ^ g (e, w) e^"" >' 2 (154) 

to operate with the evolution operator (147) on the fourier component F(u)\ , u 2 ) (see also Fig. 9, right), 
which yields 

UV(t', -oo) wa) = \c^l( Xl ) cl4(x 2 ) ^ J g r(e) 

x J^9o(e,iv) J^U(z)e- it ^-'-^e- i ^ u ' + '- e ^F(w 1 ,u^) (155) 
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and after time integrations 

U^(t',-^)F(u Jl ,oj 2 ) = \rs){rs\ J ^ e -«'(«-"- O r(e) | g 5o (e, W ) | ^ U(z) 

x w 2 ) 27T<5(w 1 — z — u) 2-k8{uoi + uo 2 — e) (156) 

Here, \ij) represents a straight (non-symmetrized) product of time-independent single-electron func- 
tions (which eliminates the factor of ±). 

If we operate on a particular energy component 

F{£)= J ^F{ui,E-ui) (157) 
the result becomes in operator form 

U^\t',-oo)F(£) = e -^ £ - H ^T(£)J^g (S,u;) 

x J^-iI{ui-u>)F(u>u£-u>i) (158) 
We can also express this result as 

U (1) (t, -oo) F(£) = e - it(£ - H °) T{£) V{£) F{£) (159) 

where 

V{£) = J^ g (£, u>) J ^ iJ(wi - w) (160) 

Here, T(£) is the resolvent (140) and go{£, w) is the operator (141). The corresponding effective in- 
teraction is obtained from the relation (68) by taking the time derivative at t = 0, which eliminates the 
resolvent, 

= pv(£)p (i6i) 

With the explicit form of the interaction (146) the matrix elements of the potential for the exchange 
of a single covariant photon becomes [64, App. A] 



(r*|V(£)|tu) = (rs I /(fc)dfc 



£ - e r - s u - (k - vy) r £ - e s - e t - (k - ij) s . 
where the A r = A sgn(e r ). The function f(k) is in the Feynman gauge given by [64, Eq. 77] 

sin(fcri 2 ) 



tu) (162) 



~4tt 2 



4:7T Z 7*12 
2 00 

= -^(l- ai . a2 )^l + ^Mkri)ji(kr 2 )C (k \l)-C (k \2) (163) 



(=0 



where j; are spherical Bessel functions and 

C (k) 

spherical tensors, closely related to the spherical 

harmonics. 

Summarizing, the contribution to the wave operator (65) from single-photon exchange, when oper- 
ating to the right on a function of the type (157), becomes 



nW(£) = Q[/«(0, -oo) = T Q (£)V(£) 



(164) 



NRC Canada 



28 

where 



T Q (£) = QT(£) 



£-H 

and the corresponding contribution to the effective Hamiltonian 



H$(£)=PV(£)P 



Can. J. Phys. Vol. , 

(165) 

(166) 



.3. Separable two-photon exchange 

(See Ref. [64, Sect.5.2.1 and App.A.2]) 

Next we consider the separable two-photon exchange for which there is an intermediate time t = t" 
with no free or uncontracted photons, as illustrated in Fig. 10 (left). The evolution operator (25) can in 
this case be expressed 

C/ (2) (t, -oo)P = U w (t,t n ) C/ (1) (f ,-oo)P (167) 

Here, the intermediate states run over all states — in the Q as well as the P space — and when the 
intermediate state lies in the model space (P), (quasi)singularities may occur. These singularities are 
removed in the reduced evolution operator (54) by including counterterms (56) 

U(t)P = U(t)P + C(t)P (168) 

We also recall the definition of the model-space contribution (MSC) (58) 

U(0)P = U(0)P + MP (169) 

where U is the evolution operator (76) with no intermediate model-space state, in this case 

U^)p = TVTqVP (170) 

The counterterm (56) is in the present case given by the product of two single-photon contributions, 




t = t' 



ft' 



- U) 



t I 

p 

F(£>) 



ft ft 



X 



■0 V 



P 

F{£) 



ft 



t = 



Fig. 10. Graphical representation of the separable two-photon-photon ladder diagram (left). This diagram is sep- 
arable, if there exists a time (represented by the dotted line) at which there is no uncontracted photon, i.e., a time 
after the first photon has been absorbed and before the second has been created. The corresponding counterterm 
(right) is a product of two operators, which evolve independently from possibly different states of the model space. 



as shown in Fig. 10 (right) 

C (2) (t) = -U {1) {t, -oo) • PU w (0, -oo) P = -U {1) {t) ■ PC/ (1) (0) (171) 
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using the notation introduced in subsection 4.2. The two factors evolve independently from (possibly 
different) states in the model space, which is indicated by the "dot". The counterterm eliminates the 
singularity, but there may be a finite remainder, which we refer to as the model-space contribution 
(MSC) (169). We shall first consider this part. 

We assume that we operate to the far right on a function of the type (157) of energy £, and that the 
intermediate model-space state has the energy £'. Using the first-order result (159), we can express the 
second-order evolution operator (167) as 

U {2) {t)P = c - it{£ - Ha) U {1) (0, 6) U {1) (0, £)P (172) 

and the counterterm — with the first factor evolving from the intermediate state — as 

C {2 \t) = -e-^ £ '- H ^ U^{0,£') ■ PU^(0,£)P (173) 

We note here that the time derivative for U as well as C eliminates the denominator of the leftmost 
resolvent. 

The MSC now becomes 

MP = (V (1) (0,£) - [/ (1) (0,£')) • PU (1) (0,£)P (174) 

Using the result (78), the last factor is 

PUW(0,£)P = V(£)P = * PV^(£)P = * PUW(£)P 

with &£ = £'- £, and with «/« = C/W (0, £') - C/« (0, £) we have 

MP= (V (1) (0,£)-?7 (1) (0,£')) •PU {1) {0,£)P= *PU (1) P (175) 

(The asterisk is used only for clarity. It notifies the position of a "fold" in the graphical representa- 
tion [47], but has no other special significance. It will mainly serve as a reminder of the position of a 
cancelled singularity, which — as we shall see — requires certain precautions.) With the definition (80) 
the MSC can be expressed 

6£ 



MP = * H$> (176) 



(Note that — and = .) The complete second-order reduced evolution operator 
(169) then becomes 



U( 2 )(0)P = U^(0)P+^-*H^ (177) 

The result above is exact also for the quasi-degenerate case. The difference ratio can be expanded 
as discussed in Appendix .5 

^£- = ^£- + 2^r S£ + v.^r S£ (178) 

which in the limit of complete degeneracy yields 



' ■• r 1 ■• l> + l *^-*II™ ( 17 9) 
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The second-order contribution to the wave operator (66) then becomes 
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n^p = qu {2 \o)p 



S£ 



cff 



d£ 



cir 



(180) 



where £1 is the wave operator (90) without intermediate model space states. 

The second-order contribution to the effective interaction is obtained by means of the relation (68). 
Since the expression (180) is valid only for t = 0, it can not be used to evaluate the time derivative. 
Instead, we have to use the original definition (168), and using the expressions (172) and (173), we find 



H$ = PV{£) Y(£) VP - PV{£') Y P {£) V{£)P = H^> + 



- m 



SH, 



eff 



5£ 



(i) 

eff 



(181) 



where = PVYqVP is the H - bar operator (80) with no intermediate model-space states. The 
last term is by definition the model-space contribution, which appears in this order only for energy- 
dependent interactions. In the case of complete degeneracy the difference ratio tends to the derivative, 
as before, 



H. 



(2) 



eff 



-L- 



8£ HeS 



U eS + 



dH. 



(i) 



(i) 



d£ cff 



(182) 



.4. Separable three-photon exchange 

The treatment of the exchange of three separable covariant photons is quite analogous to the previous 
case. From the expansion (87) we have 

fjwp = fj(3 )p + (pmpuwp _ uwp . uWp) + (pwpuwp - uwp ■ U {2) P) 

+ (U (1) PU (1) P - U {1) P • U (1) P) (PU {1) P - -PU W P) (183) 
By generalizing the result of the preceding Appendix we obtain the relation 



AP BP-AP- BP 



SA 
~5£ 



* B 



(184) 



where A is an arbitrary operator and B can be U, U or U. Using this relation, the second and third 
terms above become 



(pwpuwp-uw-puwp) 
(uWpump-uWp-uWp) 



S£ 
SUV 
8£ 



*PU (2) P 



sow (2) 
5£ eff 



(185) 
(186) 



In the last term in Eq. (183) we have to apply the rule (184) twice, yielding 

(pwpuwp - U {1) P ■ U {1) P) (PUMP - -PUMP) 

,)l ' * pu^p (puwp - .puwp) = 4 i 5 -^- * pu^p) * p&wp 

o£\ o£ ) 



S£ 

S /SUM 



= — ( '-^-^ * HQ ) * H 



r(l) 

6£ V S£ ' ^ cff J " cff 
From the previous Appendix (Eq. 178) we have 

SU^ SUM 



S£ 



S£ 



+ S£\ S£ H ° s ) 



(187) 
(188) 
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and the complete result then becomes 

£/(3)p = {/(3)p + 



S£ 



SUM 

ss 



*H 



(2) 
eff 



31 



(189) 



This is an exact expression in this order, also for a quasi-degenerate model space. In the case of com- 
plete degeneracy this becomes 



In terms of the ft operators the results above then become 



d£ 



*H. 



(2) 
off 



Q(3)p = 0(3)p + ^*^) + ^*^) 



S£ 



S£ 



off 



(190) 



(191) 



In order to obtain the third-order effective interaction, we consider the time derivative of the relation 
(183) (only the first factor is time dependent). This yields 



H 



(3) 



S<3) . ^< (1) 



' i: ^ 11% + S M*H$ + if^S- * H«) * (192) 



r l 



'Off - "off T" J£ "Off ■ 

which using the relation (182) can be expressed 



cff S£ V <5£ 



F (3) = #3) 



'off 



+ 



(i) 

off 



S£ 



<5# 



(2) 



off 



8£ 



*H 



(i) 

off 



(193) 



.5. Expansions 

We have seen above that when there are multiple singularities, it is important to take the difference 
ratios before the singularities are removed. We shall illustrate this here by a simple mathematical ex- 
ample. 

We consider a function f(x) of the variable x. We define the first-order difference ratio 



Sf_ = 5 XOiX f = f(x) - f(x ) 
Sx 



8x 



X — X 



(194) 



which can be expanded in a Taylor series 



Tx~ = %r = f'( x o)+l.r(x )(x-x )+^.r'(x )(x-x ) 2 - 



i 

4! 



f IV (x )(x-x f + --- (195) 



where 



/'(so) 



dx 



(196) 



X — Xq 



etc. 



Similarly, we define the second-order difference ratio 

1 1 



^ 2 / ^rr'rr &x ,xf _ _ r. 

S 2 x Sx Sx 

1 



x' — X 



+ ^f IV (xo) [(x' - x Q ) 2 + (x' - x Q )(x - x ) + (x- x ) 2 ] + ■ 



(197) 
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the third-order difference ratio 

3x"x' Sx'x 3x ,xf 1 ,.„/, \ 1 elV I \f i / i n o \ i 

~ / Fo) + jj/ (x )(x + x + x -3x ) + 



<5 3 x (5a; <5x Sx 
the fourth-order difference ratio 



3!' 



(198) 



(199) 



and so on. 

Generalizing these results, we have in the limit, when the differences tend to zero 



S n f 1 d n f 
S"x n\ d n x 



(200) 



This relation is frequently used in the present paper. 
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